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Robert  J.  Hinde 
University  of  Tennessee 


1.  Project  objectives  and  summary  of  accomplishments 

Specific  impulse  (/sp)  is  a  standard  performance  metric  for  rocket  propulsion  systems, 
a  metric  that  measures  the  momentum  transferred  to  the  rocket  per  unit  weight  of  propel¬ 
lant  consumed.  For  the  current  state-of-the-art  liquid  oxygen-liquid  hydrogen  propellant 
system,  7sp  =  389  s  [Sut86].  However,  if  the  hydrogen  propellant  contains  A1  atoms  at 
a  mole  fraction  of  5%,  the  theoretical  Isp  of  this  propellant  system  increases  by  9%  to 
425  s;  if  the  hydrogen  propellant  contains  5%  B  atoms,  the  theoretical  7sp  jumps  to  470  s, 
for  an  improvement  of  21%  [Car93].  These  potentially  dramatic  increases  in  performance 
sparked  considerable  interest,  within  both  the  Air  Force  Research  Laboratory  and  the  Air 
Force’s  contractor  community,  in  the  development  of  techniques  for  synthesizing  gram-scale 
quantities  of  atom-doped  solid  H2  with  dopant  mole  fractions  near  5%. 

To  guide  and  evaluate  these  synthetic  techniques,  it  is  imperative  that  researchers  have 
a  method  for  measuring  accurately  the  concentration  of  atomic  dopants  in  solid  H2  matrices 
with  thicknesses  on  the  order  of  one  to  several  millimeters.  For  atomic  dopants  with  allowed 
electronic  transitions  in  the  visible  and  ultraviolet  spectral  regions,  such  as  A1  and  B  atoms, 
electronic  spectroscopy  can  quickly  disclose  the  presence  or  absence  of  a  dopant.  However, 
these  allowed  electronic  transitions  become  optically  thick  in  millimeter-thick  H2  solids 
with  atomic  dopant  concentrations  above  roughly  0.5%  [Faj93];  consequently,  electronic 
spectroscopy  cannot  provide  quantitative  concentration  information  in  the  concentration 
range  of  interest. 

Atom-doped  solid  H2  matrices  also  absorb  weakly  in  the  infrared  (IR)  spectral  re¬ 
gion  [Faj98].  These  dopant-induced  IR  absorption  features  arise  from  nominally  IR  inac¬ 
tive  Q^O)  pure  vibrational  transitions  of  individual  H2  molecules  in  the  matrix,  transitions 
that  become  IR  active  in  the  presence  of  dopants  [Hin03].  Because  these  dopant-induced 
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features  axe  inherently  weak  [Faj98],  1-mm  thick  samples  of  doped  solid  H2  will  remain 
optically  thin  (with  transmittance  values  of  0.1  or  higher)  even  at  dopant  concentrations 
approaching  5%.  This  makes  these  features  attractive  prospects  for  measuring  dopant 
concentrations  in  thick  H2  solids. 

The  objective  of  the  research  documented  in  this  final  performance  report  was  the 
development  of  methods  for  computing,  from  first  principles,  the  integrated  intensity  of 
these  dopant-induced  Q1(0)  vibrational  transitions  in  solid  H2  matrices  containing  atomic 
substitutional  impurities.  Progress  toward  this  goal  was  made  on  several  fronts  for  the 
case  in  which  the  impurity  is  an  S-state  atom  such  as  an  alkali  metal  atom  or  a  rare  gas 
atom:  (1)  The  mechanism  by  which  S-state  atomic  impurities  generate  IR  activity  in  the 
vibrational  coordinate  of  nearby  H2  molecules  was  elucidated.  (2)  A  theoretical  framework 
was  established  for  computing  the  line  shapes  and  integrated  intensities  of  dopant-induced 
Qj(0)  features  in  solid  H2  matrices.  (3)  Path  integral  Monte  Carlo  simulations  of  rare- 
gas-doped  solid  H2  were  performed.  Based  on  these  simulations,  testable  quantitative 
predictions  were  made  for  the  dopant-induced  Q^O)  line  shapes  in  solid  H2  matrices  con¬ 
taining  Ar,  Kr,  or  Xe  atoms.  (4)  Ab  initio  calculations  for  Ar-H2  dimers  were  performed 
to  provide  input  for  future  calculations  of  the  intensity  of  Ar-induced  Qx (0)  features  in 
Ar-doped  solid  H2. 

2.  Atom-induced  infrared  absorption  features  in  solid  hydrogen 

The  central  accomplishment  during  the  performance  period  covered  by  this  final  report 
is  the  development  of  a  theoretical  framework  for  understanding  and  interpreting  the  pure 
vibrational  H2  Q1(0)  features  induced  by  S-state  atomic  substitional  impurities  in  the 
infrared  (IR)  absorption  spectra  of  solid  hydrogen  matrices.  This  theoretical  framework 
includes:  (1)  an  explanation  of  the  mechanism  by  which  atomic  impurities  generate  Qx(0) 
IR  activity  in  neighboring  H2  molecules,  and  (2)  a  method  for  relating  the  line  shape  of  the 
dopant-induced  Qx(0)  feature  to  the  vibrational  dependence  of  the  dopant-H2  potential 
energy  surface. 

Before  describing  this  theoretical  framework,  we  first  review  some  of  the  relevant  prop¬ 
erties  of  pure  solid  H2  [Sil80,  Van83].  In  solid  H2,  the  rotational  and  vibrational  motions 
of  individual  H2  molecules  are  essentially  unperturbed  by  the  interactions  between  adja¬ 
cent  molecules;  hence  the  rotational  and  vibrational  quantum  numbers  j  and  v  of  each 
molecule  remain  good  quantum  numbers  in  solid  H2  [Sil80,  Van83].  The  theoretical  frame¬ 
work  summarized  in  this  final  performance  report  assumes  that  these  quantum  numbers 
also  remain  good  in  solid  H2  matrices  that  have  been  doped  with  atomic  substitutional 
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impurities.  This  will  be  the  case  if  the  dopant-H2  interaction  is  only  weakly  anisotropic 
and  only  weakly  dependent  on  the  H2  vibrational  quantum  number. 

The  ground  state  of  the  solid  hydrogen  crystal  is  solid  parahydrogen  (pH2),  a  hexag¬ 
onal  close  packed  lattice  [Van62]  with  nearest-neighbor  spacing  of  3.79  A  [SilSO]  in  which 
each  H2  molecule  has  the  rotational  quantum  number  j  =  0.  The  rotational  degrees  of 
freedom  of  these  j  =  0  H2  molecules  are  described  by  the  isotropic  spherical  harmonic 
Y00{6,  (j>)\  the  pH2  crystal  is  thus  composed  of  spherically  symmetric  objects  with  mul¬ 
tipole  moments  that  are  identically  zero,  and  is  held  together  by  dispersion  forces.  The 
H2  molecules  in  the  pH2  crystal  can  therefore  be  treated  as  point  particles,  each  with  an 
embedded  vibrational  degree  of  freedom.  The  translational  zero  point  motion  of  individ¬ 
ual  molecules  in  the  solid  H2  matrix  is  rather  large,  so  that  the  width  of  the  first  peak  in 
the  radial  distribution  function  for  solid  H2  is  approximately  18%  of  the  nearest-neighbor 
distance,  even  in  the  T  =  0  K  limit  [Sil80]. 

A.  The  induction  mechanism 

Consider  the  Q^O)  vibrational  excitation  ( v,j )  =  (1,0)  (0,0)  of  a  single  H2 

molecule  in  solid  pH2.  Upon  excitation,  this  molecule  induces  pairwise  transition  dipole 
moments  with  each  of  its  neighbors,  as  indicated  schematically  in  Fig.  1(a);  these  transi¬ 
tion  moments  arise  from  short-range  overlap  interactions  between  pH2  molecules  [McK74, 
McK89,  Mey89].  The  high  symmetry  of  the  pH2  crystal  causes  these  pairwise  transition 
moment  vectors  to  sum  to  zero,  as  is  clear  from  Fig.  1(a).  The  Qi(0)  transition  is  therefore 
IR  inactive  in  the  pure  solid. 

However,  if  one  of  the  pH2  molecules  in  the  solid  is  replaced  with  an  atomic  sub¬ 
stitutional  impurity,  the  lattice  symmetry  is  broken  in  the  vicinity  of  the  dopant.  This 
symmetry  breaking  may  arise  both  from  physical  restructuring  of  the  pH2  crystal  lattice 
in  the  vicinity  of  the  impurity  (structural  symmetry  breaking)  and  from  the  fact  that  the 
dopant-pH2  interactions  generate  a  Q^O)  transition  moment  that  differs  from  the  one 
associated  with  pH2-pH2  pairs  (electrical  symmetry  breaking).  In  either  case,  the  net 
transition  moment  for  Qx(0)  excitations  localized  on  molecules  neighboring  the  impurity 
is  nonzero,  as  shown  in  Fig.  1(b),  and  the  vibrational  coordinates  of  these  molecules  are 
therefore  IR  active.  Atomic  impurities  thus  induce  IR  absorption  activity  in  pH2  molecules 
in  the  host  matrix  as  a  direct  consequence  of  the  intermolecular  interactions  between  the 
impurities  and  the  pH2  molecules. 

Because  the  overlap  interactions  that  give  rise  to  the  Q^O)  transition  moment  are 
short-ranged,  the  pH2  molecules  that  are  in  direct  contact  with  the  atomic  impurity  are  the 
only  ones  that  contribute  directly  to  the  intensity  of  the  Qj(0)  transition.  This  observation 
plays  a  key  role  in  our  analysis  of  the  dopant-induced  Qa(0)  feature  below. 
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B.  The  vibrational  Hamiltonian  of  atom-doped  solid  parahydrogen 

Here  we  consider  dopant-induced  Qx(0)  absorption  features  in  a  crystal  that  contains 
N  pH2  molecules  and  an  S-state  atomic  impurity.  The  position  vector  R0  denotes  the 
impurity’s  position  and  the  vectors  {Rn  :  n  =  1,2,3, . . . ,  JV}  denote  the  centers  of  mass 
of  the  H2  molecules  in  the  crystal. 

Because  the  vibrational  motion  of  the  pH2  molecules  occurs  on  a  much  faster  time 
scale  than  any  of  the  system’s  other  motions,  we  make  the  simplifying  assumption  that 
the  H2  vibrations  can  be  adiabatically  separated  from  the  other  coordinates;  we  will  use 
vk  to  denote  the  vibrational  quantum  number  of  the  H2  molecule  at  Rfe.  The  crystal  is 
therefore  described  completely  by  the  position  vectors  {R0,  Rx,  R2, . . . ,  R^  }  and  the  list 
of  vibrational  quantum  numbers  {v1,v2,  • .  • ,  vN}. 

If  the  H2-H2  and  dopant-H2  potentials  do  not  depend  on  the  H2  vibrational  coor¬ 
dinates,  then  the  translational  degrees  of  freedom  {R0,  Rx,  R2, . . . ,  R^}  decouple  from 
the  H2  vibrational  coordinates.  In  this  limiting  case,  the  wavefunction  of  the  doped 
pH2  solid  can  be  written  as  the  product  of  a  translational  wavefunction  T(Q),  where 
Q  =  {R0,R1,R2,...,Rn},  and  a  wavefunction  $  that  describes  the  vibrational  degrees 
of  freedom  of  the  N  H2  molecules.  When  the  H2-H2  and  dopant-H2  potentials  depend 
only  weakly  on  the  H2  vibrational  coordinates,  this  separation  of  translational  and  vi¬ 
brational  degrees  of  freedom  remains  a  good  approximation  to  the  doped  solid’s  true 
wavefunction.  Henceforth,  this  approximation  is  assumed  to  be  valid,  and  the  transla¬ 
tional  wavefunction  T(Q)  is  assumed  to  be  independent  of  the  list  of  vibrational  quantum 
numbers  {v1,v2,  •  •  ■  ,vN}- 

The  initial  state  |i)  of  the  doped  pH2  crystal,  before  the  IR  absorption  event  occurs, 
is  one  in  which  each  H2  molecule  in  the  system  is  in  its  vibrational  ground  state.  After 
IR  excitation,  the  crystal  is  in  a  final  state  |/)  whose  vibrational  wavefunction  $  can 
be  described  within  an  excitonic  formalism  [Fre31]  as  a  linear  combination  of  localized 
excitations  residing  on  individual  H2  molecules.  These  localized  excitations  are  denoted 
| n,  v),  which  represents  a  localized  state  in  which  the  vibrational  quantum  number  of  pH2 
molecule  n  is  vn  =  v,  while  the  vibrational  quantum  number  of  all  other  pH2  molecules  is 
zero. 

For  the  dopant-induced  Qx(0)  excitation,  the  crystal’s  final  vibrational  wavefunction 
$  can  be  written  as  a  linear  combination  of  the  v  =  1  manifold  of  localized  excitations. 
Hence  the  final  state  |/)  is 

|/>  =  *(Q)  x  ^2Cn\n,v  =  l)  .  (1) 

n 
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The  coefficients  {(?„}  that  define  the  crystal’s  final  state  following  IR  excitation  can  be 
obtained  by  diagonalizing  the  excitonic  Hamiltonian  matrix  H  describing  the  vibrational 
degrees  of  freedom  of  the  collection  of  N  pH2  molecules. 

The  Hamiltonian  operator  H  of  the  doped  pH2  crystal  depends  on  the  translational 
variables  Q  and  the  H2  vibrational  quantum  numbers  {vx,v2, . . .  ,vN]  through  the  voca¬ 
tionally  averaged  potential  energy  surface  V(Q;v1,u2, . . .  ,vN)  for  the  solid.  The  matrix 
elements  of  this  operator  in  the  basis  of  localized  vibrational  excitations  can  be  written 


Hn'n  — 


(n>=1l[/ 


1)U2)  ' 


,V)|*(Q)|2dQ 


(2) 


and  are  seen  to  contain  integrals  over  the  probability  density  function  |T(Q)|2  for  the 
translational  zero  point  motion  of  the  doped  solid. 

The  diagonal  Hamiltonian  matrix  element  Hn  n  is  simply  the  energy  of  the  localized 
excited  state  |n,t>  =  l)  in  which  the  H2  vibrational  excitation  resides  on  molecule  n.  The 
off-diagonal  elements  Hn,  n  are  averages  over  ^(Q)]2  of  the  terms  in  the  potential  energy 
surface  V(Q;  iq,  v2, . . .  ,vN)  which  couple  the  vibrational  coordinates  of  pH2  molecules  n 
and  n'. 

The  eigenvectors  of  H  describe  the  final  states  accessed  in  the  IR  absorption  event, 
as  indicated  in  Eq.  (1).  The  coefficients  {Cn}  that  describe  a  particular  final  state  |/)  are 
the  elements  of  a  column  eigenvector  of  H.  The  absorption  intensity  for  this  final  state  is 
proportional  to  the  square  of  the  corresponding  transition  dipole  moment  integral,  which 
for  dopant-induced  Qx(0)  excitations  is 


<*|M|/>  =  / M„(Q)|*(Q)|2  dQ  . 

n  ** 


(3) 


In  Eq.  (3),  M„(Q)  is  the  Q1(0)  transition  moment  associated  with  excitation  from  the 
initial  state  |z)  to  the  localized  final  state  | n,v  =  l).  This  transition  moment  depends, 
through  Q,  on  the  structure  of  the  doped  pH2  crystal  in  the  neighborhood  of  the  dopant. 
Finally,  the  dopant-induced  Q1(0)  IR  absorption  line  shape  for  a  doped  pH2  solid  can  be 
obtained  by  weighting  each  final  state  |/)  by  its  absorption  intensity  (i|M|/)2. 

The  integrals  in  Eqs.  (2)  and  (3)  can  be  evaluated  using  quantum  many-body  tech¬ 
niques  such  as  diffusion  quantum  Monte  Carlo  (DQMC)  and  path  integral  Monte  Carlo 
(PIMC)  simulations  [And75,  And99,  Cep96,  Cha81,  Suh91].  These  techniques  provide 
computationally  efficient  methods  for  treating  the  large-amplitude  correlated  motion  of 
neighboring  molecules  in  the  doped  pH2  solid.  An  accurate  treatment  of  this  large- 
amplitude  motion  is  important  for  a  number  of  reasons.  First,  the  transition  dipole  mo¬ 
ments  M  in  Eq.  (3)  are  dominated  by  overlap  interactions  between  the  dopant  and  its  pH2 
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neighbors.  These  transition  moments  depend  strongly  on  the  dopant-pH2  distance  [Mey89] 
and  become  substantially  more  important  when  they  are  averaged  over  the  underlying  zero 
point  motion  [HinOO].  Second,  as  Fig.  1  indicates,  the  transition  dipole  that  controls  the 
intensity  of  the  dopant-induced  IR  absorption  features  is  a  vector  quantity  that  depends 
simultaneously  on  the  positions  of  the  pH2  molecule  undergoing  excitation,  its  impurity 
neighbor,  and  its  eleven  pH2  neighbors.  A  faithful  accounting  of  the  correlated  motions  of 
these  thirteen  molecules  is  thus  critical  for  computation  of  the  transition  moment. 

C.  A  tight-binding  model  for  the  vibrational  Hamiltonian 

A  simplified  model  for  the  vibrational  Hamiltonian  matrix  H  of  doped  pH2  solids  can 
be  developed  by  extending  the  tight-binding  vibron  Hamiltonian  introduced  to  explain  the 
IR  and  Raman  line  shapes  of  pH2  solids  containing  j  —  1  H2  molecules  as  dopants  [Fel95, 
Sea64,  Van59].  In  this  tight-binding  vibron  Hamiltonian,  the  vibrational  coordinates  of 
nearest-neighbor  (nn)  pH2  molecules  are  assumed  to  be  coupled  through  the  H2-H2  in- 
termolecular  potential;  pH2  molecules  which  are  not  nearest  neighbors  are  assumed  to  be 
remain  decoupled.  Under  the  assumption  of  only  nearest-neighbor  coupling,  the  tight- 
binding  vibron  Hamiltonian  operator  can  be  written  as  [Fel95] 

H  =  ^TEk\k,v  =  l)(fc,u  =  l|  +  A  \k,v  =  l)(n,v  =  l|  (4) 

k  fc,n=nn 

where  |n,  v  =  l)  again  represents  the  state  of  the  crystal  in  which  the  vibrational  excitation 
is  localized  on  pH2  molecule  n.  In  this  equation,  Ek  is  the  Q^O)  transition  energy  for  the 
localized  vibrational  excitation  associated  with  pH2  molecule  k\  this  transition  energy  can 
be  obtained  from  Eq.  (2)  by  setting  n'  —  n  =  k.  Intermolecular  interactions  in  the  solid 
reduce  Ek  from  its  gas  phase  value  of  Ek  =  4161.1  cm-1  to  Ek  =  4152.2  cm-1  [Van68]. 
The  parameter  A  in  Eq.  (4)  represents  the  vibrational  coupling  between  nearest-neighbor 
molecules.  Analysis  of  the  IR  absorption  spectrum  of  pairs  of  j  =  1  H2  molecules  isolated 
in  solid  pH2  [Zha98]  indicates  that  the  hopping  parameter  A  =  -0.225  cm-1  for  a  nearest- 
neighbor  pair  of  H2  molecules. 

In  a  doped  pH2  crystal,  the  Q1(0)  transition  energy  Ek  for  pH2  molecules  next  to 
an  impurity  will  differ  slightly  from  that  for  pH2  molecules  distant  from  the  impurity 
because  of  differences  between  the  H2  vibrational  dependence  of  the  H2-H2  and  dopant-H2 
potentials.  We  account  for  this  by  introducing  a  parameter  A E  that  quantifies  the  shift  in 
the  Q1(0)  transition  energy  for  pH2  molecules  adjacent  to  the  dopant.  For  these  molecules, 
Ek  =  4152.2  cm-1  -  A E\  for  all  other  pH2  molecules,  Ek  remains  at  4152.2  cm-1.  This 
leads  to  a  one-parameter  family  of  vibron  Hamiltonians  defined  by  the  parameter  A E. 
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The  matrix  elements  expressed  in  Eq.  (2)  are  then  equivalent  to  matrix  elements  of  Eq.  (4) 
evaluated  in  in  a  basis  set  consisting  of  the  localized  vibrational  excitations  \n,v  —  l). 
Diagonalizing  this  matrix  and  then  weighting  the  eigenvalues  by  their  respective  absorption 
intensities  ultimately  yields  the  dopant-induced  (^(0)  line  shape  as  described  in  section  B 
above. 

Figure  2  shows  the  normalized  line  shapes  for  several  values  of  A E,  ranging  from 
0.25  cm-1  to  2  cm-1;  the  calculations  leading  to  the  results  shown  in  Fig.  2  are  described 
in  detail  in  [Hin03].  At  low  values  of  the  detuning  parameter  A E,  the  dopant-induced 
Q^O)  absorption  feature  shows  substantial  IR  intensity  over  the  entire  pH2  vibron  band 
ranging  from  v  —  4149.5  cm-1  to  v  =  4153.1  cm-1.  As  A E  increases,  the  dopant-induced 
absorption  feature  shifts  to  the  red  and  sharpens;  for  A E  =  2  cm-1  the  feature’s  full  width 
at  half  maximum  is  only  0.02  cm-1.  This  evolution  in  the  shape  of  the  dopant-induced 
feature  as  the  detuning  parameter  increases  reflects  a  transition  from  a  delocalized  IR- 
active  vibron  at  small  A E  to  a  strongly  localized  IR-active  vibron  at  large  A E.  When 
Aj E  is  small,  the  IR-active  pH2  molecules  adjacent  to  the  dopant  can  readily  share  their 
excitation  with  neighboring  pH2  molecules;  consequently,  the  vibrons  that  contribute  to 
the  IR  absorption  spectrum  can  be  quite  spatially  delocalized.  As  A E  increases,  however, 
it  becomes  more  and  more  difficult  for  a  vibrational  excitation  localized  in  the  dopant’s 
first  “solvation  shell”  to  hop  away  from  the  dopant,  and  the  vibrational  coordinates  of 
pH2  molecules  in  this  solvation  shell,  which  are  the  only  IR-active  molecules  in  the  solid, 
become  effectively  decoupled  from  the  rest  of  the  system. 

3.  Detuning  parameters  for  rare-gas-doped  solid  hydrogen 

If  the  detuning  parameter  A E  for  a  particular  atomic  dopant  is  known,  the  vibron 
Hamiltonian  defined  by  Eq.  (4)  allows  us  to  predict  the  line  shape  of  the  dopant-induced 
Qj  (0)  absorption  feature  in  pH2  solids  containing  this  atomic  dopant.  We  have  used  PIMC 
simulations  to  estimate  the  A E  values  for  Ar,  Kr,  and  Xe  dopants  in  order  to  provide  line 
shape  predictions  for  these  dopants  that  can  eventually  be  compared  with  experiment. 

The  A E  value  for  molecule  k  in  a  rare-gas-doped  pH2  solid  can  be  estimated  as 

A  £(fe)  =  J  AF(fc)(Q)|^(Q)|2  dQ  (5) 

where  the  quantity 

A7(fc)(Q)  =  V(Q;  wj  =  0,  u2  =  0, . . . ,  vk_,  =  0,vk  =  0,  vk+1  =  0, . . . ,  vN  =  0) 

-V(Q;uj  =0,u2  =0,...,vfc_1  =0,ufe  =  l,vk+1  =0,...,vN  =  0) 
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is  the  negative  of  the  change  in  the  potential  energy  of  the  doped  pH2  solid  when  molecule 
k  makes  a  transition  from  its  v  =  0  vibrational  level  to  its  v  =  1  vibrational  level.  (Recall 
that  A E  is  defined  above  to  be  a  positive  number  if  the  Q2  (0)  transition  of  a  H2  molecule 
red  shifts  under  the  influence  of  a  neighboring  dopant.)  In  principle  this  function  can  be 
decomposed  into  contributions  from  H2-H2  and  dopant-H2  potentials;  however,  at  the 
moment  there  are  no  accurate  H2-H2  potentials  that  explicitly  incorporate  vibrational 
dependence  of  the  individual  H2  molecules.  Consequently,  we  approximate  AVW(Q)  as 
follows.  We  first  assume  that  the  8.9-cm-1  red  shift  of  the  H2  Q^O)  transition  in  pure 
solid  pH2  can  be  partitioned  among  an  individual  molecule’s  12  nearest  neighbors  in  an 
additive  fashion.  If  molecule  k  is  in  the  first  solvation  shell  of  the  rare  gas  dopant,  it 
has  11  pH2  neighbors  and  is  therefore  missing  (1/12)  x  8.9  =  0.74  cm-1  of  this  red  shift. 
Consequently  we  define 

AK^(Q)  =  I^op— H2 v  =  0)  —  Kiop— H2(^0’ v  =  -0  ~  0*74  cm  (7) 

where  Vrdop_H2(R0,  Rk]  v)  is  the  rovibrationally  averaged  dopant-pH2  potential  for  a  pH2 
molecule  in  vibrational  level  v.  Recall  that  the  vectors  R0  and  Rfc  specify  the  locations  of 
the  dopant  and  pH2  molecule  k,  respectively. 

The  integral  in  Eq.  (5)  is  computed  for  doped  pH2  solids  at  T  =  4  K  using  discrete¬ 
time  PIMC  simulations  [Cha81]  that  employ  the  primitive  Trotter  factorization  of  the 
Boltzmann  factor  exp (-H/kT).  We  model  the  rare-gas-doped  pH2  solid  as  a  collection 
of  179  pH2  molecules  and  one  dopant  atom  at  a  number  density  equal  to  that  of  pure 
solid  pH2  at  T  =  4  K;  the  particles  are  initially  placed  on  a  hexagonal  close  packed  lattice 
containing  six  layers  of  30  particles  each  [Che96].  We  use  periodic  boundary  conditions  to 
eliminate  surface  effects.  The  potential  energy  of  the  doped  solid  is  taken  to  be  a  pairwise 
additive  sum  of  H2-H2  potentials  [Sil78]  and  dopant-H2  potentials  [LeR87].  We  assume 
that  the  pH2  molecules  obey  Boltzmann  statistics. 

Figure  3  shows  how  the  A E  value  for  pH2  molecules  surrounding  an  Ar  dopant  atom 
depends  on  the  Trotter  number  NT.  We  see  that  A E  converges  to  the  value  A E  = 
0.20  cm-1  for  Trotter  numbers  iVT  >  48.  For  our  production  calculations,  we  use  the 
Trotter  number  iV-p  =  64.  Our  final  estimates  of  the  A E  values  for  Ar-,  Kr-,  and  Xe- 
doped  pH2  are:  A E  =  0.20  cm-1  for  Ar  dopants,  A E  =  0.85  cm-1  for  Kr  dopants,  and 
A£  =  1.68  cm-1  for  Xe  dopants.  The  statistical  uncertainties  in  these  estimates  are 
roughly  ±0.02  cm-1. 

These  detuning  parameters,  when  interpreted  in  light  of  the  results  summarized  in 
Fig.  2,  indicate  that  Ar-doped  solid  pH2  should  exhibit  a  broad,  symmetric,  featureless 
dopant-induced  Qx(0)  spectrum  spanning  the  entire  pH2  vibron  band;  this  is  in  good 


8 


agreement  with  experimental  observations  [Hin02].  Kr-doped  solid  pH2  should  exhibit  an 
asymmetrical  dopant-induced  Q^O)  feature  with  peak  intensity  near  v  =  4150.5  cm-1, 
while  Xe-doped  solid  pH2  should  exhibit  a  fairly  narrow  dopant-induced  Q1(0)  feature 
near  u  =  4149.75  cm— 1,  with  a  low-intensity  tail  to  the  blue. 

From  our  PIMC  simulations,  we  can  extract  probability  distribution  functions  for  the 
distance  between  the  dopant  atom  and  its  nearest-neighbor  pH2  molecules.  These  are 
shown  in  Figure  4.  Each  panel  in  Fig.  4  also  includes,  for  reference  purposes,  the  nearest- 
neighbor  distance  distribution  for  pure  solid  pH2  at  T  =  4  K;  this  is  computed  using  PIMC 
methods  as  well. 

4.  Ab  initio  calculations  of  the  Ar-H2  dipole  moment 

The  integrated  intensities  of  dopant-induced  Q^O)  IR  absorption  features  in  atom- 
doped  solid  pH2  are  determined  by  the  transition  dipole  moment  integral  defined  in  Eq.  (3). 
This  integral  in  turn  depends  on  the  Qx(0)  transition  moment  Mn(Q)  associated  with 
excitation  from  the  initial  state  |i)  to  the  localized  final  state  |n,u  =  l).  If  the  dipole 
moment  operator  for  the  doped  pH2  solid  can  be  written  as  a  pairwise  sum  of  interaction- 
induced  dipoles  (see  Fig.  1),  then 

M„(Q)  =  £m<‘,“)(R*'RJ  (8) 

k 

where  is  the  dipole  moment  generated  by  interactions  between  particles  k  and  n, 

and  the  sum  in  Eq.  (8)  is  over  only  the  twelve  nearest  neighbors  of  pH2  molecule  n. 

If  n  is  a  pH2  molecule  and  k  is  either  another  pH2  molecule  or  an  S-state  atomic 
dopant,  then  k  and  n  are  both  spherically  symmetric  particles.  The  dipole  moment  function 
then  depends  only  on  the  interparticle  vector  R^.n  =  Rn  —  R*..  Furthermore,  in 
this  case  rrdfc,n^  must  point  along  the  interparticle  vector  Rfc  n,  so  that  we  have 

=  R  (9) 

where  mSk'n\Rk  n)  is  now  a  scalar  function  of  the  interparticle  distance  Rk<n-  The  de¬ 
pendence  of  this  function  on  the  vibrational  coordinate  of  pH2  molecule  n  contains  the 
essential  information  needed  to  predict  the  integrated  intensities  of  dopant- induced  Qj  (0) 
features  in  doped  pH2  solids. 

In  the  case  that  k  and  n  are  both  pH2  molecules,  the  vibrational  dependence  of  this 
function  has  been  computed  using  ab  initio  quantum  chemical  methods  [Mey89].  Similar 
calculations  are  needed  for  dopant-pH2  dimers  to  enable  evaluation  of  the  transition  dipole 
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moment  matrix  elements  defined  by  Eq.  (3).  We  have  used  ab  initio  quantum  chemical 
techniques  to  compute  the  interaction-induced  dipole  moment  M  for  the  Ar-H2  dimer. 
This  dipole  moment  is  the  derivative^  of  the  Ar-H2  interaction  potential  V  with  respect 
to  the  magnitude  of  an  externally  applied  homogeneous  electric  field  F:  M  =  (dV/dF)F=0. 
The  interaction  potential  V  is  computed  using  Gaussian  [Gau03]  at  the  CCSD(T)  level  of 
theory  employing  the  standard  counterpoise  correction  [Boy70]. 

Our  frozen-core  calculations  employ  standard  aug-cc-pVTZ  atom-centered  basis  sets 
that  have  been  extended  with  an  additional  set  of  diffuse  functions.  These  additional 
functions  consist  of  one  primitive  Gaussian  of  each  symmetry  with  an  exponent  one-third 
as  large  as  the  smallest  exponent  of  that  symmetry  in  the  standard  aug-cc-pVTZ  basis 
sets.  Midway  between  the  Ar  nucleus  and  the  H2  center  of  mass,  we  place  a  (3s3p2d)  set  of 
bond  functions  [Tao92]  to  improve  the  description  of  dispersion  interactions.  Only  valence 
electrons  are  correlated  in  our  ab  initio  calculations.  Extensive  tests  show  that  core-valence 
correlation  makes  an  insignificant  contribution  to  the  interaction-induced  dipole  moment 
of  the  Ar-H2  dimer. 

We  describe  the  configuration  of  the  Ar-H2  dimer  using  Jacobi  coordinates,  in  which 
R  is  the  Ar-H2  distance,  r  is  the  H-H  bond  length,  and  7  is  the  angle  between  the  Ar-H2 
vector  and  the  H-H  bond.  For  specificity,  we  place  the  dimer  in  the  (x,  z)  plane  with  the 
Ar  atom  at  the  origin  and  the  H2  center  of  mass  on  the  positive  z  axis;  when  0  <  7  <  7r/2, 
the  H  atom  which  is  closer  to  Ar  has  a  positive  x  coordinate. 

The  derivative  ( dV/dF)F=0  is  estimated  from  a  finite- difference  analysis  of  the  in¬ 
teraction  energies  obtained  when  the  dimer  is  perturbed  by  external  fields  of  strength 
F  =  ±0.001  au.  Such  an  analysis  gives  accurate  dipole  moments  if  the  interaction  energy 
can  be  expressed  as  a  quadratic  function  of  the  external  field  strength  for  small  values  of 
F;  tests  show  that  this  is  indeed  the  case  for  field  strengths  less  than  or  equal  to  0.003  au. 
Hence  the  finite-difference  analysis  employed  here  should  give  reliable  dipole  moments. 

For  the  configuration  ( R  =  6.856  a0,r  =  1.401  a0,7  =  0)>  we  obtain  a  field-free 
frozen-core  Ar-H2  interaction  energy  of  V  =  -52.50  cm-1;  this  is  in  excellent  agreement 
with  the  value  of  -52.56  cm-1  obtained  from  a  frozen-core  CCSD(T)  calculation  [Woo98] 
employing  the  aug-cc-pV5Z  basis  set  for  H  and  the  d-aug-cc-pV5Z  basis  set  for  Ar.  This 
demonstrates  that  bond  functions  can  effectively  compensate  for  inadequacies  of  atom- 
centered  Gaussian  basis  sets  in  computing  the  binding  energy  of  weakly  bound  dimers. 

(i)  In  Gaussian  03  revision  B.05,  the  orientation  of  an  externally  applied  field  F  is  defined 
so  that  when  a  molecule’s  permanent  dipole  moment  vector  M  points  in  the  same  direction 
as  F,  the  field-dipole  interaction  raises  the  molecule’s  energy. 
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Furthermore,  we  find  that  the  interaction-induced  dipole  moment  M  changes  by  less  than 
1%  when  our  extended  triple-zeta  basis  sets  are  replaced  by  either  larger  aug-cc-pVQZ 
basis  sets  or  smaller  aug-cc-pVTZ  basis  sets.  This  indicates  that  the  basis  sets  used  here 
are  large  enough  to  yield  converged  interaction-induced  dipole  moments  at  the  CCSD(T) 
level  of  theory. 

At  a  fixed  Ar-H2  distance  R,  the  Cartesian  components  of  the  interaction-induced 
dipole  moment  can  be  written  as  [Mey86a] 

oo  oo 

MX  =  C2j,fcP2j(cos7)(r  -  r0)k  (10) 

jfc=0  j= 1 

and 

M*  =  £  £  D2j,kP2j(COSl)(r  -  ro)k  (n) 

*=0  j= 0 

where  Pn  is  a  Legendre  polynomial,  P™  is  an  associated  Legendre  polynomial,  and  r0  is  a 
reference  H-H  bond  length  which  we  choose  to  be  1.4  a0  here.  Figure  5  shows  the  angular 
dependence  of  Mx  and  Mz  for  Ar-H2  at  R  =  5  a0  and  r  =  1.4  a0.  The  ab  initio  results 
are  accurately  reproduced  by  a  fit  that  includes  the  j  =  1  term  in  Eq.  (10)  and  the  j  =  0 
through  j  =  2  terms  in  Eq.  (11).  To  determine  the  number  of  powers  of  (r  -  r0)  that  must 
be  retained  in  Eqs.  (10)  and  (11),  we  calculate  Mx  and  Mz  for  Ar-H2  as  functions  of  r  at 
R-  5  a0  and  three  values  of  7;  our  results  are  shown  in  Fig.  6.  We  see  that  for  r  values 
between  r  =  0.9  a0  and  2.1  a0,  truncating  these  power  series  expansions  at  (r  -  r0)2  is 
adequate. 

Figures  5  and  6  show  that  we  can  determine  the  coefficients  of  the  leading  terms  in 
Eqs.  (10)  and  (11)  for  a  given  R  value  by  computing  the  Ar-H2  dipole  moment  at  nine 
configurations,  specified  by  three  choices  for  r  and  three  choices  for  7.  Here  we  choose 
r  =  1.1,  1.4,  and  1.7  a0  and  7  =  0, 7r/4,  and  7t/2.  We  calculate  the  interaction-induced  Ar- 
H2  dipole  moment  at  these  nine  configurations  for  several  R  values  between  5  a0  and  10  a0. 
At  a  fixed  value  of  R,  we  then  integrate  Eqs.  (10)  and  (11)  over  specific  H2  rovibrational 
wavefunctions  to  obtain  rovibrationally  averaged  dipole  moments  and  rotationally  averaged 
transition  dipole  moments.  In  Fig.  7,  we  present  the  rovibrationally  averaged  Ar-H2  dipole 
moment 

(M)  oo  =  (v  =  0  ,j  =  0|  Mz\v  =  0  J  =  0)  (12) 

and  the  Ar-H2  Q^O)  v  =  1  0  transition  dipole  moment 

(M)01  =  (v  =  0,  j  =  0| Mz\v  =  1  ,j  =  0)  .  (13) 
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In  Eq.  (13),  the  relative  phases  of  the  H2  v  =  0  and  v  =  1  vibrational  wavefunctions  are 
chosen  so  that  the  integral  (v  =  0,  j  =  0|r|u  =  1,  j  =  0)  >  0. 

The  (M) oo  and  (M)01  functions  can  be  fit  reasonably  well  by  an  expression  of  the 

form 

(M)  =  D7/R7  -  c0  exp(c1i?  +  c2R 2) 

using  these  parameters  (all  in  atomic  units): 

(M)00  :  D7  =  278,  c0  =  29.815,  cx  =  -1.1693,  c2  =  -0.0543 

(M)01  :  D7  =  166,  c0  =  4.450,  c2  =  -0.9479,  c2  =  -0.0496 

Our  (M)00  curve  for  Ar-H2  is  in  qualitative  agreement  with  that  presented  in  [Mey86b], 
where  the  vibrationally  averaged  Ar-H2  dipole  moment  (M)00  was  estimated  by  perform¬ 
ing  ab  initio  calculations  for  Ar-H2  dimers  with  the  H2  fragment  frozen  at  its  v  =  0 
vibrationally  averaged  bond  length.  If  we  evaluate  the  rotationally  averaged  (j  =  0)  Ar- 
H2  dipole  moment  at  this  H-H  bond  length  using  Eqs.  (10)  and  (11)  and  our  computed 
C  and  D  coefficients,  our  results  are  roughly  10%  larger  than  those  presented  in  [Mey86b] 
for  Ar-H2  distances  R  between  5  and  6  a0 . 

5.  Publications 

The  following  publications  report  work  supported  in  part  or  in  full  by  AFOSR  grant 
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•  “High-resolution  spectroscopy  of  HC1  and  DC1  isolated  in  solid  parahydrogen:  direct, 
induced,  and  cooperative  infrared  transitions  in  a  molecular  quantum  solid”,  D.  T. 
Anderson,  R.  J.  Hinde,  S.  Tam,  and  M.  E.  Fajardo,  J.  Chem.  Phys.  116,  594-607 
(2002). 

•  “Probing  quantum  solvation  with  infrared  spectroscopy:  infrared  activity  induced  in 
solid  parahydrogen  by  N2  and  Ar  dopants”,  R.  J.  Hinde,  D.  T.  Anderson,  S.  Tam, 
and  M.  E.  Fajardo,  Chem.  Phys.  Lett.  356,  355—360  (2002). 

•  “Infrared-active  vibron  bands  associated  with  substitutional  impurities  in  solid  para¬ 
hydrogen”,  R.  J.  Hinde,  J.  Chem.  Phys.  119,  6-9  (2003). 

•  “Mg-He  and  Ca-He  van  der  Waals  interactions:  approaching  the  Born-Oppenheimer 
limit”,  R.  J.  Hinde,  J.  Phys.  B  36,  3119-3128  (2003). 
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Figure  1.  Symmetry  breaking  induces  IR  activity  in  solid  pH2  matrices  containing  atomic 
substitutional  impurities.  Open  single  circles  represent  v  =  0  pH2  molecules;  double  circles 
represent  pH2  molecules  undergoing  a  v  =  1  «-  0  excitation,  (a)  In  a  close  packed  pure 
pH2  solid,  excitation  of  a  single  pH2  molecule  generates  twelve  transition  dipole  moment 
vectors,  six  of  which  are  shown  here.  The  twelve  vectors  sum  to  zero  because  of  the 
underlying  symmetry  of  the  close  packed  lattice,  (b)  In  a  pH2  solid  containing  an  atomic 
substitutional  impurity  (shaded  circle),  excitation  of  a  pH2  molecule  next  to  the  impurity 
yields  a  nonzero  net  transition  dipole  moment  because  the  transition  moment  for  dopant- 
pH2  dimers  differs  from  that  for  pH2-pH2  dimers. 
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Figure  2.  Theoretical  dopant-induced  pH2  Qj  (0)  absorption  spectra  for  several  values  of 
the  detuning  parameter  AE.  From  bottom  to  top,  AF?  —  0.25  cm  ,  0.5  cm  ,  0.75  cm  , 
1  cm-1,  1.5  cm-1,  and  2  cm-1.  The  spectra  are  offset  vertically  for  clarity.  Note  that  the 
spectra  for  A E  >  1  cm-1  have  been  reduced  in  intensity  for  ease  of  presentation. 
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Figure  4.  Distribution  of  dopant-H2  nearest-neighbor  distances  in  (a)  Ar-doped,  (b)  Kr- 
doped,  and  (c)  Xe-doped  solid  pH2  at  T  =  4  K.  The  dashed  line  in  each  panel  represents 
the  distribution  of  nearest-neighbor  H2-H2  distances  in  pure  solid  pH2  at  T  =  4  K. 
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Figure  5.  Dependence  of  the  Ar-H2  dipole  moment  on  7  at  R  =  5  a0  and  r  =  1.4  a0. 
Open  circles  indicate  values  of  Mx;  open  squares  indicate  values  of  Mz.  The  lines  are  fits 
to  truncated  Legendre  polynomial  expansions  as  described  in  the  body  of  the  report. 
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Figure  6.  Dependence  of  the  Ar-H2  dipole  moment  on  r  at  R  =  5  a0.  Panel  (a)  shows 
Mx  as  a  function  of  r  at  7  =  7t/4;  panel  (b)  shows  Mz  as  a  function  of  r  at  three  7  values. 
The  solid  lines  are  quadratic  fits  to  the  ab  initio  calculations. 
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